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SUMMARY 

The  purpose  of  tMs'paper  Is  to  illustrate  the  application 
of  a  result  in  matrix  theory  to  the  problem  of  determining  the 
linear  differential  equation  whose  solutions  are  the  products 
of  the  solutions  of  two  given  linear  differential  equations. 
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0N  THE  LINEAR  DIFFERENTIAL  EQUATION  WHOSE 
SOLUTIONS  ARE  THE  PRODUCTS  OF  SOLUTIONS 
OF  TWO  GIVEN  LINEAR  DIFFERENTIAL  EQUATIONS 

By 

Richard  Bellman 


§1 .  Introduction 


It  mbe  observed  by  Newton  that  a  simple  way  to  obtain 
the  power  series  expansion  for  the  function 

( 1 )  u  -  (arc  sin  x  ) 2 

was  to  form  the  second  order  linear  differential  equation 

(2)  ( 1  — x 2 ) ^ ^  -  x  —  -  2, 

dx"  dx 

and  eolve  It  by  means  of  a  power  series  expansion. 


In  another  direction,  It  was  shown  by  Appell,  [l]  ,  that 
If  u, ,  u0  represent  two  linearly  Independent  solutions  of 

(>)  u  r  p(t)u  ♦  q(t)u  =  0, 

2  2 

tnen  u^  ,  u  u0,  uQ  represent  three  linearly  Independent  solu¬ 
tions  of  the  third  order  linear  differential  equation 


(4)  u 
of ,  also 


Op(t)u'  '^(2p2(t  )+p'  (t)4-4q(t))u'  +  (4p(t)q(t)+2q'(t))u*0 
Whittaker  and  Watson,  l2j  . 


This  result  Is  useful  In  connection  with  the  determination 


of  the  power  series  expansion  for  the  square  of  the  hypergeometrlc 
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functior.  f(a,b,c;x)  ,  and  plays  a  role  In  the  study  of  the 
Mathleu  function,  [2]. 

There  are  several  ways  of  determining  the  equation  In 
(4),  elnce  the  problem  Is  analogous  to  that  of  finding  the 
polynomial  whose  roots  are  various  symmetric  combinations  of 
the  roots  of  a  given  polynomial. 

In  this  paper,  we  shall  present  a  new  method  based  upon 
an  Interesting  result  concerning  matrix  differential  equations. 

§  - •  Preliminary  Lemma 

The  crux  of  the  method  Is  the  well-known 
Lemma .  Let  Y  and  Z  be  respectively  the  so  .utjons  of 

(.)  j-f  -  A(t)Y,  Y(0)  -  I, 

-  ZB(t),  2(0)  -  J. 

Then  tne  solution  of 

(2)  3^  -  A(t)X  ♦  XB(t) ,  X(0)  -  C, 

Is  given  by 

(3)  X  -  YCZ . 

Verification  provides  an  Immediate  proof.  Let  us  assume 
tnat  A ( t )  and  B(t)  satisfy  the  condition  of  being  lntegrable 
over  any  finite  Interval  . 

•Whlie  this  paper  was  In  the  process  of  being  typed,  a  paper  on 
this  theme  appeared-  H.  P.  Sandr.&m,  A  Square  and  a  Product  of 
Hypergeometrl  c  Functions,  Quart.  Jour,  cl'  Math.,  Vol.  7(1956), 
PI  1  P'3- A  5^  • 


§3 •  Application 
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Let  ue  now  apply  this  resu.t  to  the  problem  of  finding 

the  3rd  order  .Inear  differentia  equation  whose  solutions  are 

2  2 

Ui'  ulu2'  u2'  "^ere  u  and  u2  are  tw0  linearly  Independent 
so lutlons  of  (1.3). 


Without  ioss  of  generality,  let  and  u?  be  determined 
by  the  boundary  conditions. 


(1) 

ux (0)  -  1  , 

uj(0)  -  0, 

u2(0)  -  0, 

u2(0)  -  i . 

1 

Setting  u  »  v,  we 

see  that  ( 

(2) 

1 

U  -  V 

1 

V  -  -p  ( t )  V  - 

q(t)u . 

Let 

(3) 


A  ( t ) 


0  i 

~q ( t )  -p ( t )  . 


Then  the  matrix  BOiutlon  of 
(<♦)  u’  -  A  ( t )  U ,  U(0)  -  I, 

le  given  by 

(5.)  U  - 


U,(t) 


and  the  solution  of 


(6) 


VA  ( t ) T  , 


V(0) 


by  v«lr ,  the  transpose  of  U 
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Prom  the  lemma  In  §2,  we  deduce  that  the  solution  of 
(7)  x'  -  A  ( t )  X  ♦  XA(t  )T,  X(0)  -  C, 

1 s  given  by 

(tt)  X  -  u  c  u  T. 


Taking  C  to  be  a  symmetric  matrix, 


(9) 


C  ■ 


1 


cA 


■  i 


we  see  that  X  is  given  by 

(10)  X  ■ 


c1u12+2c2uiu2+c3u22  ciulul+c2^u2ul^u2ui ^C^U2U2 

i  1 2  i  «  \  o 

Vl  *2c2UlVC?U2 


VciuiVc2<u2u.+Vji)+cJu3u3 


Writing 
(11)  X  - 

the  equation  In  (7)  Is  equivalent  to  the  ayttem 


fX  1  *2 

V*2  Xj  ’ 


(■2) 


2(a 


111 


u 


X2  "  a2-Xl  *  ^al.  *  a22  ^  X2  f 

X-^  "  ^  ®  p  '  x  2  +  a  °  2  x  3  ^  " 


Eliminating  x  and  x^,  we  obtain  a  third  order  linear 

differential  equation  for  x .  wnose  genera,  solution  Is 

2  2 

c  u  +  2c,u  4-  c,u0  ,  where  c  ,  c„,  and  c,  are  arbitrary 

1  .  c  ;  c  l  c  ; 


constants  . 
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Similarly,  eliminating  x^  and  x^,  we  obtain  the  equation 

n  it  ip 

wnose  general  solution  Is  c  ^  +  c-^u^  ;  eliminating 

x^  and  x^,  we  obtain  an  equation  whose  solutions  are  the  deri¬ 
vatives  of  the  solutions  of  the  equation  obtained  by  the  ellml— 

i 

nation  of  x2  and  x^. 

§4 .  The  General  Case 

In  stating  the  lemma  In  §2,  we  Ignored  any  discussion  of 
the  dimensions  of  Y  and  Z.  It  Is  clear  that  the  result  Is 
valid  If  A(t)  and  Y  are  m  x  m  matrlcee,  B(t)  and  Z  n  x  n  matrices, 
and  C  and  X  m  x  n  matrices. 

Using  the  technique  sketched  In  §3,  we  can  obtain  the 
linear  differential  equation  of  order  m  n  whose  solutions  are 
the  products  of  the  solutions  of  a  linear  differential  equation 
of  order  m  and  one  of  order  n. 
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